Ballistic side jump motion of electrons and holes in semiconductor quantum wells 
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We investigate the ballistic motion of electrons and holes in III-V semiconductor quantum wells 
with spin-orbit coupling and a homogeneous in-plane electric field. As a result of a non-perturbative 
treatment of both of these influences, particle wave packets undergo a pronounced side jump per- 
pendicular to the field direction. For wave packets of sufficient width the amplitude of this motion 
can be estimated analytically and increases with decreasing field strength. We discuss the scaling 
behavior of the effect and some if its experimental implications. 



Effects of spin-orbit coupling in semiconductors have 
become a particularly lively direction of solid-state re- 
search. Most of the activities here are embedded in the 
emerging field of spintronics. An early paradigm of this 
area is the spin field-effect transistor proposed by Datta 
and Das over fifteen years ago [1]. In this device pro- 
posal, electrons are injected from a spin-polarized source 
electrode into a quantum well in which the dominant con- 
tribution to spin-orbit interaction is given by the Rashba 
coupling, 



(1) 



Here p is the momentum of the electron confined in a 
two-dimensional geometry, and a the vector of Pauli ma- 
trices. The Rashba coefficient a is tunable by an electric 
gate across the well and can therefore be varied experi- 
mentally [2]. As a result, in the course of the electron's 
ballistic motion through the quantum well its spin un- 
dergoes a controlled rotation which can be detected by 
another spin-polarized drain electrode. To enable the in- 
jection of electrons into this device, a certain voltage has 
to be applied between source and drain electrode, lead- 
ing to a possibly small but nevertheless finite in-plane 
electric field in the quantum well. In the present let- 
ter we investigate the effect of such an electric field in 
combination with spin-orbit coupling. As we shall see 
below, the interplay between spin-orbit coupling and an 
in-plane electric field adds an important feature to the 
ballistic electron dynamics: The electron performs a side- 
jump motion perpendicular to the direction of the electric 
field. The field here is assumed to be homogeneous, in 
contrast to side jump motion upon scattering on impu- 
rities which is discussed intensively in circumstances of 
the anomalous Hall effect [3] . 

Moreover, as it was pointed out recently, spin-orbit 
coupling of the type (1) does not only lead to a rotation 
of the spin, but has also an effect on the orbital degree 
of freedom: the electron performs a zitterbewegung, i.e. 
an oscillatory motion superimposed on the nonrelativis- 
tic dynamics of the particle wave packet [4]. In the ab- 
sence of an electric field the single-particle Hamiltonian 
reads 7i = ^ -|- T-Lr where m is the effective band mass. 



In this case a fully analytical description of the zitter- 
bewegung can be given in terms of the time-dependent 
position operator r7/(i) — exp{iHt/'h)f{0) exp{—i'Ht/Ti) 
in the Heisenberg picture. As we shall see, the side jump 
phenomena to be discussed below are accompanied by 
zitterbewegung of this type. 

Another important recent development regarding spin 
dynamics in semiconductors is the prediction of the in- 
trinsic spin Hall effect [5-8] . This effect is again a conse- 
quence of spin-orbit coupling and amounts in a spin cur- 
rent (as opposed to a charge current) driven by a perpen- 
dicular electric field. Therefore, this situation is similar 
to the Datta-Das transistor with the electric field being 
provided by the source-drain voltage. Here spin orbit 
coupling of the Rashba type has the peculiarity that the 
spin Hall conductivity vanishes in an infinite system in 
the presence of any spin relaxation mechanism [9] . How- 
ever, this conclusion holds only in the thermodynamic 
limit, and for a device of mesoscopic size and sufficient 
carrier mobility, spin Hall transport signaled by spin ac- 
cumulation at the sample edges should be considered as 
possible. 

Let us consider an electron in an n-doped quantum 
well being subject to Rashba spin-orbit coupling and a 
constant in-plan electric force F, i.e. the single-particle 
Hamiltonian reads — + T^R. ^ F'^- Thus, the 
time-dependent momentum operator is given by ph (t) = 
p(0) + Ft, while the spin operators fulfill the equations 



-pUtWi)] 



(2) 
(3) 
(4) 



with (T//(0) — a. Note that the position operator r does 
not occur in the above equations, and therefore the time- 
dependent spin operators (7 nit) can be formulated as 
a function of the zero-time Schrodinger operators p, a 
only. In particular, when the operators 'Jfj{t) are ap- 
plied on eigenstates of the momentum p = p(0), the 
quantities pnit) = p(0) + Ft become real numbers, ren- 
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dering Eqs. (2)- (4) as an system of ordinary differential 
equations for the 2 x 2-matrices <JH{t)- Unfortunately an 
explicit analytical solution of these equations for general 
directions of the force F and the initial momentum p{Q) 
does not seem to be possible. However, such a solution 
can be given if F and p(0) are collinear. Since the Rashba 
Hamiltonian (1) is invariant under rotations of spin and 
momentum in the xy-plane, we can choose, without loss 
of generality, F and p(0) to point along the a;-direction 
where we find 

aUt) = <y'' cos (^-^[p^t+\Fey^ 

+ a^sin{^^{pH+\Fey^ , (5) 
<(^)=<7^ (6) 
a|,(t) = -a-sin(^^(^p-t+iFt2^) 

+ a'cos(^^(^pH + \Ft'^^ . (7) 

Note that the argument of the trigonometric functions is 
the integral dt'p^{t'), generalizing the situation with- 
out an external field and therefore constant momentum. 

The time-dependent position operators fnit) fulfill the 
equations 

|x.W = ^ + J<W, (8) 
at m n 

-T:yH[t) = — -^(^hW- (9) 

at m n 

Let us now consider the expectation values {fH{t)) for an 
initial state with momentum p\\F along the .T-dircction, 
spin in positive z-direction and (rH(0)) = 0. Here we 
have {xnit)) = [p^t + ^Ft'^)/m, and the y-component 
reads 

{VHit)) = -l[dt'sin (p^f + ^Ft'^)) . (10) 

If the initial spin direction is reversed, {yH{t)) changes 
sign. The case of vanishing initial momentum p(0) = is 
particularly interesting. Here the physical problem con- 
tains only two length scales which are conveniently cho- 
sen as the Rashba length = Yi' jmo? and the field- 
dependent length If = sJajF. The latter quantity de- 
termines the amplitude of the side jump motion to be dis- 
cussed now. Figure 1 shows a plot of {yait)) for p(0) = 0, 
F = 1.0meV//im, and a band mass m = 0.023mo, cor- 
responding to the conduction band mass on InAs, where 
mo is the bare electron mass. The Rashba parameter 
has been fixed to a = O.OleVnm which is a realistic 
value for InAs [8]. The main panel shows shows the a 
real-space plot of {fait)) with the time t ranging from 
zero to t = 50ps. One can clearly distinguish two phases 



of the electron motion: In a first phase with {xnit)) ^ 
200nm, the transverse expectation value (y//(i)) starts at 
(yj?(0)) = and undergoes a pronounced monotonic side 
jump motion to \ {yH{t)) \ ~ 60nm. In a following second 
phase {ynit)) performs a zitterbewegung, i.e. an oscilla- 
tory motion around this value with its amplitude decreas- 
ing with increasing longitudinal distance {xH{t))- The 
period of the zitterbewegung as a function of (xnit)) is 
approximately given by wXr « 1040nm [4]. The inset of 
figure 1 shows the same data for {ynit)) but as a function 
of time t. After a time interval oi t = 50ps the electron 
wave number is k^{t) = p''{t)/h = Ft/h < COSnm-^ 
which is still close to the F-point in III-V semiconductors. 
Thus, the assumptions underlying the effective-mass ap- 
proximation and the Rashba Hamiltonian (1) are per- 
fectly valid. 

As seen from figure 1, the magnitude of the side jump 
can be obtained by performing the formal limit 

lim^iVHit)} = -lj^ dtsin (^^Ft^) (11) 
= -^./|, (12) 

where we have used the integral dt sin(i^) = -^/tt/S ~ 
0.62. For the above parameters we have Ip = lOOnm 
corresponding to a side jump of about 62nm, in perfect 
agreement with the numerical plot. Note that according 
to Eq. (10) the band mass m does not enter the function 
{yH{t)), but is inversely proportional to {xH{t))- Thus, 
changing the effective mass as m m' amounts in a 
rescaling of the abscissa in the main panel of figure 1 by 
a factor of m/m'. Thus the above considerations relat- 
ing to InAs arc trivially extended to other materials, say 
GaAs. In the latter case the Rashba coupling is typi- 
cally somewhat smaller which can be compensated by a 
smaller in-plane field F, resulting the the same length 
scale If, i-e. the same magnitude of side jump. 

In the above discussion we have concentrated on plane- 
wave eigcnstatcs of the momentum. However, a localized 
electron in a quantum well is not described realistically 
by a single plane wave but rather a superposition of them. 
We therefore consider a Gaussian wave packet, 

m -^^j d^ke-'-^'^^-^^^r ( J ) , (13) 

where the initial spin direction is again along the z-axis. 
We numerically solve for the dynamics of such a wave 
packet by superimposing plane-wave solutions whose ini- 
tial momenta have been chosen from a grid centered 
around A;o in reciprocal space. Because the in-plane field 
F is in general not collinear with such initial momenta, 
these plane-wave solutions have also to be generated nu- 
merically. Moreover, in such a procedure necessarily re- 
quires a wave vector cutoff. In all numerical simulations 
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to be presented below we have carefully checked that the 
data is well-converged with respect to grid spacing and 
wave vector cutoff. 

For simplicity, let us concentrate again on the case 
of zero initial group wave number, ko = 0. Then the 
width d of the wave packet is the only new length en- 
tering the problem. Figure 2 shows numerical simula- 
tions of the real space dynamics for different wave packet 
width d with otherwise the same parameters as before. 
If = lOOnm. As seen from the figure, the dynamics of a 
single plane wave discussed in Eqs. (10), (12) is a good 
approximation to the wave-packet motion if d » Ip- 
In fact, for d = lOOOnm = IQlp the dynamics of the 
wave packet is very close to the plane-wave result shown 
in figure 1, which continues to be a good approxima- 
tion up to (i w 500nm = 51 p, whereas for smaller d 
the nature of the motion clearly changes. The criterion 
d ^ Ip follows from the observation that Eqs. (2)- (4) 
are invariant under the rcscaling F qF, p 
t '—^ t/y/q with q being a real parameter. Therefore, 
in order to ensure invariancc, the width d of the wave 
packet has to scale as rf/Vf which uniquely determines 
Ip = \JalF as the relevant length scale. Note also that 
Eqs. (2)- (4) are also invariant under the replacements 
(p*^(t),(T^(t)) ^ (-p*^(t),-a},(i)), i e {x,y}. As a con- 
sequence, for the above situation of a wave packet sym- 
metric around its group wave vector ko = and the force 
F pointing along the x-direction, we have the scaling be- 
havior {ynit)) 1-^ iUHit)) / whereas {xH{t)) remains 
unchanged. The latter conclusions hold for any symmet- 
ric wave packet, independently of its form and width; 
in particular, the magnitude of the side jump scales al- 
ways the same way as Ip. Analogous findings hold for 
the rcscaling a i-^ qa, p i-^ p/^/<li t ^/^/Q leading 
to {xH{t)) ^ {xH{t))/^, ivHit)) ^ ^{yH{t)). There- 
fore, figure 2 is of universal character since data for other 
values of a, can be obtained by a simple rcscaling of 
the coordinate axes. 

So far we have concentrated on the Rashba Hamilto- 
nian (1) as an effective contribution to spin orbit cou- 
pling. The situation is changed only marginally if the 
Rashba coupling is replaced with the Dresselhaus term 
[10] in its linear approximation [11]. When both terms 
are present an analytical solution for dnit) does not 
appear to be possible even for initial conditions with 
p{0) II F. An exception is the case when the Dresselhaus 
and the Rashba term are equal in magnitude; Here side 
jump and zitterbewegung are absent because of an addi- 
tional conserved quantity arising at this high-symmetry 
point [12]. 

Let us now turn to p-doped quantum wells of III-V 

semiconductors. Here the low-energy physics is domi- 
nated by the heavy holes, and the spin-orbit interaction 
term analogous to the Rashba coupling on electron spins 
reads [13] 



^R = i^{P-^+-pW) , (14) 

using the notations p± = Px i ipy, a± = ± ia^ . The 
Pauli matrices operate on the total angular momentum 
states with spin projection ±3/2 along the growth direc- 
tion chosen as then z-axis,and a is the spin-orbit coupling 
coefficient. This Hamiltonian contains two length scales 
which are, similarly to the previous case of conduction- 
band electrons, aptly chosen to be the field-dependent 
length Ip = {a/FY/^ and the Rashba length = 
mo? jfi' where m is the effective heavy- hole band mass. 
Moreover, an analytical plane-wave solution can again be 
found if the initial momentum p(0) is coUinear with the 
in-plane force F . Choosing this direction, again without 
loss of generality, to lie along the x-axis, one finds 

<7^(t) = cos £ dt' (p^ + Ft'f^ 

^cT^sinl^^ jyt'ip^+Ft'fy (15) 
a^(^)=a^ (16) 
<7|f (0 = a- sin (^^ 1^ dt' (p- + Ftf'^ 

+ cos (^^ dt' {p^ + Ft'f^ , (17) 

where the remaining integral is elementary. For a plane- 
wave state with p{0) = 0, (rf/ (0)) = one finds again 
{xnit)) = Ft^/2m and 

iVHit)) = -P dt'{Ft'f sin {^,FH'') . (18) 
The magnitude of the side jump is given by the limit 

^lim (yH(t)) = -3 • 2^/4 j j dxx'^ sin (a;*) 

« -1.428- [p.. (19) 

Thus, the length scale Ip plays a similar role as Ip 
for conduction-band electrons. Moreover, analogously 
to this previous case, the behavior of plane-wave states 
is mimicked by proper wave packets provided d ^ Ip, 
a conclusion following from similar scaling arguments 
as before. The situation is illustrated in figure 3 with 
numerical simulations for values for m, a appropriate 
for heavy holes in p-doped GaAs quantum wells and 
Ip = lOnm. For wave packet width up to d « 5lp 
the magnitude of the side jump is quantitatively well de- 
scribed by Eq. (19), whereas for smaller values of d the 
nature of the dynamics changes. 

In summary, the pronounced side jump found in the 
dynamics of electrons and holes in semiconductor quan- 
tum wells is the non-perturbative effect of the inter- 
play of spin-orbit coupling and a homogeneous electric 
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field. The magnitude of the side jump motion is gov- 
erned by the field-dependent length scales Ip = \fojF 
and If = {a/F)^^'^ for electrons and holes, respectively, 
provided that the width d of the particle wave packet 
being initially at rest fulfills d ^ 5{If,If)- A possible 
scenario for the experimental detection of the side jump 
effect are two opposite quantum point contacts of width d 
whose centers are displaced by A 0.2d. The conductiv- 
ity of this setup should be maximal if the experimentally 
tunable parameters a (a) and F are adjusted such that 
hilF) ~ A. 

Another interesting question is the connection between 
the ballistic side jump motion and the intrinsic spin Hall 
effect in diffusive systems. In a heuristic but appealing 
picture one can interpret spin hall transport as iterated 
side jumps occurring during the ballistic motion between 
scattering events. In fact a theoretical synthesis of ex- 
trinsic and intrinsic mechanisms of spin Hall transport 
[8] was put forward in Rcf. [14] where the side jump con- 
tribution to the spin Hall conductivity was found to be 
independent of disorder and particle interactions. 

I thank I. Adagidcli and D. Weiss for useful discus- 
sions. This work was supported by the SFB 689 "Spin 
Phenomena in reduced Dimensions" . 
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FIG. 1. Main panel: real-space motion (rW(t)) for a plane 
wave with p(0) = 0. The inset shows {yH{t)) as a function of 
time t. 
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FIG. 2. Real-space motion {rH{t)) for electron wave pack- 
ets with zero initial group velocity and different width d. 
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FIG. 3. Real-space motion (rff(t)) for hole wave packets 
with zero initial group velocity and different width d. 
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